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INTRODUCTION 


T h i s w a r k d e a 1 s w i t h a t h e o r e t i c a 1 s t u d y o f t !" > 

r a d i a t i a n / a n n e a 1 :i. n g r e s p c :< n s s o f M 0 S3 I C s „ A ]. t h a u g hi m a n y 

experiments have been per -farmed in this field, no comprehensive 
theory dealing with r ad i at i on /anneal i ng response has; been 
proposed . Many attempts have been made to apply linear response 
theory in several works 1 ” 4 , but no theoretical foundation has been 
pres; anted 

T h e 1 i n e a r r e s p o n s e t h e o r y o u 1 1 i n e d h e r e i s c a p a b 1 e o f d e s c r i b i n g 
a b r o a d a r e a o f r a d i a t i o n / a n n e a I .i n g r e s p o n s e p h e n o rn e n a i n M !I) S I C s , 
i n p a r t i c u 1 a r , b o t h s i m u 1 1 a n e o u s i r r a d i a t. i. a n a n d a n n e a 1 i n g , a s 

w e 1 1 a s s h o r t. - a n d 1 o n g t e r m a n n e a 1 i n g , i n c 1 u d i n g t h e c a s e w h e n 

anneal ing is nearing completion . For the first time, a simple 

p I'" o c e d are i s d e v i s e d to d e termi n e t he r e s p o n s e f u n c t i on from 

experimental r ad i at i on /anneal i ng data. In addition, this 

p r o c e d u r e e n a b 1 e s u s t o s t u d y t h e e f f e c t s o f v a r i a b 1 e t e mpera t u r e 
and dose rate, effects which are of interest in spaceflight. 

In the past, the shift in threshold potential due to 
r a d i a t i o n / a n n e a 1 ;i. n g * ~ + h a s u s u a 1 1 y to e e n a s s u m e d t o d e p e n d o n o n e 
variables the time lapse between an impulse dose and the time of 
observation (Eq. (1))., While such a suggestion of uniformity in 
1 1 m e i s c e rtainly tru e f o r a b r o a d r a n g e o f r a d i a t ion / a n n e a 1 in g 
p h e n o rn e n a , i t m a y n o t h o 1 d f o r s o m e r a n g e s o f t h e v a r i a b les o f 
interest (temperature, dose rate, etc. ) , We therefore propose 

here a response function which is dependent on two variables: the 

time of observation and the time of the impulse dose. This 
dependence on two variables allows us to extend our theory to the 
t r e a t m e n t o f a v a r i a b 1 e d o s e r a t e . 

Final ly, the 1 inear theory is generalised to the case in which the 
response is nonl inear with impulse dose, but is proportional to 
some i mpul se function of dose. A method to determine both the 
impulse and response functions is presented. 
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I „ GENERAL THEORY 


We begin with the assumption that the incremental response of a 
ii 0 S d e v i c e t o r a d i a t i o n / a n n e a I i n g , t h a t i s , t h e s h ;i. f t o f t h e 
t h r e s h o I d p o t e n t i a 1 , 8 V , :i. s 1 :i. n e a r w i t In r e s p e c t t o a n i m p u 1 s e 

radiation dose, 6'D, and with respect, to some function which 
depends on the time elapsed between the time of irradiation, t, 
and the time of observation, t ’ , (Fig. 1) : 

SK> ( t , t '' ) = R ( t ■■■" '!: ) SD < t ’ ) , ( 1 ) 

Eq. ( 1 ;) may be rewritten as 

6V £ t. , t ’ ) = R ( t t ■" ) D ( t ■ ) St < 1 a ) 

where R ( t-t ’ ) is the response function,, t is the time of 
observation , t.* is the time at which the radiation impulse occurs, 
S‘D is the radiation dose impulse and D"dD/dt, D being the dose 
rate . 

The response function, R, in Eq. (1) is assumed to depend on the 
t. i m e 1 a p s e b e t w e e n 1. r r a d i a t i. o n a n d o b s e r v a t i o n ,, :i. n o t h e r w o r d s , R 

i s a f Li n c t i o n o f o n e v a r i a b 1 e . T h I s i s c e r t a i n 1 y v a 1 i d f o r t h e 
case of constant dose rate,, We will also assume that this is 
v a 3. i d f a r t h e c: a s e o f a s 1 o w 3. y v a r y i n g d o s e r a t e . W e w :i. 1 1 r e f e r 
to this as the "adiabatic" case. Since the case of varying dose 
rates has not been studied so far, we will defer for the? time 
b e i n g a n e a c t d e f i n i t i o n o f w h a t. i s m e a n t by "si o w 1 y v a r y i n g 1 ’ 
with respect. to dose rate., For a very rapidly varying dose rate 
Eq. (1) obviously does not hold true and should be expressed by: 



r t‘+U' * 


F i g u re 1 . I m p u 1 s e d a s e . 

To find an expression for the total shift of the threshold 
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potential, i5'V. 
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OF POOR QUALITY, 


we integrate over the 


time of i r r a d :i. a t i o n , 


6-y. tot (ti - j d v < t , tr >.. 


( 1 c: ) 


Substituting Eq. (.1) in differential, form, we obtain 


JVtct(t) = J R (t-t’ ) dD < t ) „ 


( 2 ) 


where t is not necessarily equal to t r - „ 
t r , E q . ( 3 ) y i e 1 d s 


I n p a r t i c : u 1 a r , f o r 1:. 


J'Vt o-t. ( t ) = j R < t - 1 ’ ) d I) ( t ) 


Der ben wick and Sander 1 proposed an integral for the case of t=tr-, 
which is only valid for the case of simultaneous irradiation and 
annealing. The validity of Eq. (2) extends not only to this case,, 
but also the case in which t /t,.- , which represents pure annealing 
after the end of the irradiation. In addition,, removing this 
limitation <t~t r -) allows us to gain insight into the response 
■function,, R(t), and to find out how it can be determined. 


T h e R e s p o n s e F' u n c t i o n 


To gain insight into the meaning of the response function, R(t 
we need to study the case in which t>> tv-,, Since, in Eq„ (2), t 
t.-,, then t : '<< t and R (t--t. ) a ; R(t)„ Eq„ (2) then becomes 

t„... 


Wtot ( t ) a; J R ( t ) dD C t -' > = R ( t ) D tB t , ( 3 

0 

where the total integrated dose within the interval from t 0 to 
t ( F ig. 2) is 

Dtot = D(t r .) - I)(0) = D ( t r . ) „ 

F o r t. h e c ase of v a r i a b 1 e do s e rate, D to t is; 

t K- 


bf. rat J D ; t ■' ) d t ■' „ 



'tf n J u‘uM i^-MhSn 


If D is cons l ant , then D fe o* Dt 

D 


i 

F i g . 2 . Tot a 1 dose. 

'Thus, Eq. (3) gives us the needed insight into the? meaning of the 
r e s p onse function!; n a m e 1 y , t h a t R ( t. > is the shift o f t h e t. h r e s h o Id 
potential per unit dose? when t>> t „ In other words, to 
determine R(t) for any time, t, all we need is to measure the 
s h i f t o f t h e t h r e s h old p c:s t e n t i a 1 , ( t ) . Then, 

R(t) ■ JV to t(t)/D to t(t r ) f (4) 

pi r o v i d e d t h a t t > > t r - . 

T h u s , t h e 1 i n e a r response of M 0 S d e v i c: e s is t o t a 1 1 y d e s c r i b e d b y 
Eq. (2) . Of course, we have to know an analytical or, at least, a 
numerical expression for the linear response function R<t.)„ The 

phenomenol ogi cal .linear response theory of Eq. (1) and (2) does not 
provide us with the response function; however Eqs. (2) and (4) can 
serve as a useful tool to plan appropriate experiments and to 
c a n s t r Li c t the response f u n c t i o n u s i n g e x p erimen t a 1 d a t a . "!" o f in d 

the response function, we only need to measure the response to the 
im p u i s e a n d substitute it i n t o E q . (4) . 

A valuable technique which we shall use below in detail is the 
treatment of the case of small time, t .. In this case, the 

response function may be expanded in a power series of t (see 

Eq. (7)). Substituting this into the general Eq„ (.2) and 

integrating, we obtain the response, £ V ( t ) , as as power series of 
t .. I h e n , c o m p a r i n g t h i s expans i o n w i t h experim e n t a 1 d a t a , w e c a n 

d e t e r m i n e t h e e x p a n s ion coe-f f i c i e n t s , i , e . , t h e v a 1 u e o f R ( t ) a n d 

its derivatives (first, second, etc.) at t = 0. 

This will be discussed later. We now need to comment on the 
d e f i n i t :i. o n o f s mall time, t „ W e a r e d e a 1 i n g h e r e w i t. h a k i n e t i c: 
process and it can be characterized by a characteristic time, t ra 
(or more than one of such parameters) . Thus, the condition, 
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C-B = In ft) - In ft. /t ra ) = ].n ft,,) 


t DC > 


b e d e t e r n j i n e d 


It is clear from Fig,, 4 that the constant A cai 
easily, since it is the slope and the constant, B , is the y- 
intercept. B is actually arbitrary, since it depends on the units 
chosen for time measurement . t.„ may be a. function of many 

variables (manufacturing process, temperature, dose rate, etc,. ) 
and is therefore difficult to estimate in advance, and the value 
of C depends on the value of t„. Thus, both C and i„ remain 
unknown,, We will present, below a way to determine the constants 
t , C a n d A f r o m t h e s x p e r i m e n t a 1 d a t a ., 


A close look at Lq„ (5b) reveals that it is valid only for large t, 
.i. . e ,, , t //to. Indeed, as t approaches zero, In ft) approaches 
infinity,, which is a physically meaningless situation.. The 
simplest way to extend Eq. (5b) to small t is to substitute 1 + 
t/t,, for t/t„ in the argument of the logarithm function:; 

R(t) = A In ( 1 + t/t.,,) - c, (6) 

wh i 1 e ,, f or t > >t„ , we st i 1 1 h a ve 


R ( t ) = A 1 n ( t / 1 „ ) - c „ 


We have s o far obtained Eq. (!5b) as an expression for 
basis of examination of the experimental data,, This 
the condition of long-term annealing (t>>t r - )„ On the 
this expression was understood to be valid only for t> 


RCt) on the 
was based on 
other hand , 

4 - 

c::» •« 


1 1 w a s t h e n 1 o g i c a 1 t a g e n e r a 1 i z e t. hi i s e q u a t i o n 
the region of small t, i.e., t comparable? to or 
still much greater than t r - . 


y e x t e n d i n g i t. t o 
less than t. ,, , but. 


Eq. (6) will be directly compared to the 
flV/D can be equated to R(t) (Eq. (6) ) , 
the magnitude of t is comparable to that 
p 1 a t 1 e d w i t h r e s p e c t. t o 1 n ( 1 + 1 ) , r a t h e r 


experimental data, .:. 
u n d e r t w o c a n d .1 1 i on s s 
of t„ and 2) the data 
than to 1 n ( t ) . 


1 ) 
a r e 


L' u t o ft Fat: t o r 

A 1 1. h o u g h E q „ ( 5 b ) cl e s c r i b e s p u r e 
toiit and t > >t r ,, it can be seen 
t h e 3. o g a r i. t h m f u n c: t i o n , a n d 

wi thout 1 i mi t , 


annealing for large t in the sense 
t h a t , a s t c: o n t i n u e s t o i n c r e a s e s , 
thus R(t), continues to increase 


The unre 

al i st i c ohy 

15 

i cal s: 

i. tuat i on 

of 

a con 

tinuous i 

n crease 

i n 

anneal i ng 

wi thout 1 i m 

d. 

t. can 

be cor i' 1 

•ect 

ed by 

:i. ntroduci 

ng a. cut 

of f 

factor, f 

(t/t m ) , such 


that 








•f 

< f /f ) 

' (.. .• L. ^ l 

j. *f Q }•"' 

t<t. 

m and 






•f 

/ 4.. / ,L. V 

/ c. m , 

approach 

105 

(**) f"‘i K" 

■b > 1- 

V. .... %... m „ 



This new 

constant , 


t , i s 

def i ried 

l_ Q 

be at 

1 east an 

order 

o f 

magrti tude 

1 arger than 


1 0 „ 
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There are 
'simplest c 


many possible forms for We offer 

ut of f f unc t i on : 


her 


t h e 


1 / !I 1 + ( t / 1 } 2 3 


6 a > 


Th e 

?n 

the mod if 

i. ed 

reap 

onse 

f Ul 





R ( t. 

) = 

IIA 1 

n < t 

We 


requi r 

e t. 

hat 

t m > 

>t 0 

sue! 

1 0 .< 

;t 

< t 

.... <... rr» «( 

i. n w 

h i c h 

Eq. 

( 6 b ) 

t. 

ver 

• v 

1 arge 

t, 

we 

have 

Eq. 

(6b 

i. n 

!*! 

h i c hi , 

f c^r' 

very 

1 ar 

ge t 

, tl 


(6b) 


a region such that. 
: E q .. (6) .. Thus, for 


This satisfies the condition that. R <t) be a reversible process,. 
We shall treat the case of an irreversible process elsewhere. 

A way to determine a value for t. m is to equate the response in 
E q . ( 6 ) t o z er o i . e . , 


A 1 n ( t / 1. ) - C 


U , 


F r o m t h i s w e o b t a i n 


In 
da' 


I- = +• pC/i a 

m v. tr M 

the next section, we will determine C and A from 
and will f ind the values of both t., n and t c „ 



( 6 c ) 

ex per ;i. mental 


Simultaneous irradiation and annealing, tr-=t . 

Substituting Eq. (6) far the response function into Eq. <2a) , we 
a b t a i n t h e following expression for the s h i f t. o f t h e t h r e s h o 1 d 
potent i al : 

£V(t) = Dt-CA(H-to/t> ln(l+t/to)-(A+C) D. (7) 

To simplify the study of the time dependence of the response, it 
is convenient to substitute the shift per unit total dose for the 
shift in threshold potential: 


SV (t> /Di 


A ( 1 + 1 „ / 1 ) 1 n ( 1 + 1 / 1 0 ) - ( A + C ) 


( a ■ 


where 


Dt 
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Bot 

h 

Eqs „ 

( 7 ) 

an 

d ( 

B ) c a n b e s 

i mp .1 

i f i ed 

i n 

the c a 

se of 1 a 

rge 

t i me 

fi 

t > > 

■t 

F 

or the 

s h i 

ft per unit. 

d o s 

e, we 

ob t 

a i n 2 











< t ) /Dtot m 

A In 

<t/t„) 

.... 

(A+C) . 



( 8 a 

) 

Thi 

v~:. 

ex pr 

E’SSi o 

n 

is 

very con ven 

i ent 

t o cl e 

cl ]. 

wi th „ 

It. is s 

i. mi 1. 

a r t 

a 

Eq. 

(5 

b ) 

f or 1 

on 

9“t 

erm anneal i 

n g . 

We rnu 

st 

r erne mb 

sr, hows 

0 K~ 

tha 

t 

Eq- 

(5 

a) 

g i ve 

s 

the 

total shif 

t i n 

three 

hoi 

d pot 

ent i al 

wh 

ere a 

t::: 

Eq. 

(8 

a) 

gi ves 

th 

e shift per 

un :i. 

t tot a 

] (jj 

0 1::> 0 » 

For th 

i s 

shi f 

t 

< t > 

>t 

) , 

we ob 

ta 

i n ; 
















5V ( t ) = Dt 

A .!. n 

<t/t c> ) 

< a+c ) :i . 



( /a 

) 

It 

X s 

ear™ 

;y to 

f i 

nd 

that ,5V (t) 

reaches a 

mi n i mum 

at t-t,„, 

wh i c h :i. 

3 

the 

cn 

o 1 ut 

. i on a 

if 

the 

f ol 1 owi ng 

6? qua 

t :i. on s 








'I !-, /+ ’I 


and 




v ) 


( 1 0 ) 


It. would be interesting to investigate whether it is possible to 
reach this minimum experimentally bet ore the process becomes 
i rr ever si b 1 e .. 


P u r e a n n e a 1 i n g , t > t r - „ 

We now substitute Eq. <6) into Eq. (2) . It is obvious that 
t_ t t t_ 


J R (t-t ■" ) dD <t’ ) = . 
0 0 


R ( t-t ' ) cJD ( t " ) 


R ( t t ; ‘ ) dint’ > . 


i ) 


Making use of this, we obtains 

<5V(t> = F ( t > - F ( t—tr- ) , 
w h e r e F ( t. ) i s t. h e r i g h t -■ h a n d s i d e o f E q .. < 7 ) . 


1 :l. ) 


II. CALCULATION OF THE RESPONSE FUNCTION FROM EXPERIMENTAL DATA. 

As we have already mentioned above, the parameters A, C and t. c;> , 
which totally define the response function in Eq. (6), should be 
determined from the experimental data. These parameters are 
expected to depend on the conditions of irradiation, e.g. , the 
temperature and dose rate. They may also be affected by the 
process parameters during manufacture, particularly the oxide 
anneal temperature. but also possibly by dopant species and 
concentration and other parameters In short, the quality of MOS 
devices is "reflected" in these constants. 
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1 1 i s n o t e n o u g h si m p 1 y t o d e t e r m i n e a n d cat a 1 o g t hese c o n s t. a n t s 

f o r v a r i o u s M 0 S d e v i c 0 s , b u t r a t h e r t o s t u d y e p e r i m e n t a 1 1 y t h e 

dependence of these parameters on both temperature and dose rate 
and compare the values obtained with the parameters at the 

m a n u lac t u r i n g p r o ces s . W i t h t h is i. n m i n d , we c o m e t. o t h e 

realization that there is a strong need -for a fast, simple way to 
d e t. erm i n e t. h e s e cons t ants. 


We propose a method based on a study of the linear response for 


t<<to, 

when 

R ( t > 

can be 

e p a n d e d i n a p 0 w e r 

series of t: 


R ( t. ) 

= R ( 0 

) + R ’ 

( 0 ) t. + l/2R"(0)t a + 1 

/3R-"(0)t 3 + . .. (12) 

where 

R’ (0) 

and 

R" (0) 

ar e , respect i vel y , 

t. h e first and s e c 0 n d 


derivatives of R with respect to t. At the point t-0, R(0) is 

R ( t.=0) „ 

In particular, for t h e In ( t ) d e p e n d e n c: e , 

R ( t ) fs •- C + A ( t / 1 o ) - 1 / 2 A ( t. / 1 » , ( 1 3 ) 

or , comp ar i n g E q . ( 1 3 ) an d E q . ( 1 2 ) , 

R ( 0 ) - -C, R* (0) = A/to, R” (0) = -A/t 0 2 and R : ' " (0) = 2A/t„®. (14) 

W e n e x t s u b s t i t u t e E q . ( 1 2 ) i n t. o E q < 2 ) o r ( 2 a ) , d e p e n d i n g o n 
whether we are dealing with the case of pure annealing after 
irra d i a t i o n o r s i m ultaneo u s i r r a d i a t. i o n an d a n n e a .ling,. Aft e r 
integration, this yields an expression for the shift in threshold 
potential per unit dose given in the following forms 

6'V (t ) /Dtot = a 0 + ait + a at. 52 + a 3 t 3 . .. (15) 

l he coefficients, a tl , ai, a.-.?, etc. are determined from the 

e p e r i m e n t. a 1 data. F r om t; h e s e e perim e n t a 1 1 y deter in i n e d v a .1 u e s , 

w e w i 1 1 d e l: e r m i ne the c o n s t a n t s A , C and t 


E per - i m e n t. a 1 P rocedure 

We used data obtained from experiments performed at the GSF-C 
R ad i at ion Faci 1 i t. y and c ompared th e s e r esul ts t o o u r t h e o r y „ I n 
these experiments, RCA 10 6 rad -hard CMOS ICs were irradiated with 
6 °Co gamma rays with a dose rate equal to approximately 130 
rads /min. An IC was irradiated -for several hours, taking 
measurements of the threshold potential (evaluated at. a drain 

current, of 300 jjA) for one n channel and one p channel transistor 

e v e r y ten min u t e s „ 


Sim u 1 1. a n e o us Irradiatio n a n d Ann e a 1 in g 
S u b s t i. t. u t i n g E q ( 1 2 ) i n t. o E q „ ( 3 ) , 


we obtain 
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t. 


6V ( t ) = DR ( 0 ) J 

(J 

t; 


t 


dt. ’ + DR ( 0 ) J ( t t J ) dt : ' + 

0 

t 


+ 1 / 2 DR" ( 0 )J 

0 


(t-t’ )*dt’ + 1 / 3 DR ? ( 0 ) J ( t -t ) 3 d t p + 

0 


Integrating, we obtain 
cW(t) = DtotCR(O) + 1 / 2 R’ ( 0 )t - 1 / 6 R" ( 0 ) t 5 *} 
where D to t = D(t) = Dt . 


J4R " CO) t :3 


W e obtain t h e f o 1 1 o w i n g e p r essi o n f o r t he s h i f t of t h r e 
p ot en t. i a I b y sub st i t ut i n g E q . ( 1 4 ) i n t. o E' q . (16): 

£V(t) = D tot [-C + l/2A(t/t 0 ) -- 1/6A (t/t,,,) = - l/i2A(t/t c> )®3, 

where t<<t«». 

It is more convenient to deal with the shift per unit dose in 
■form of Eq. (15), where the constants are 

a — C , a i := A / 2 1 0 , a s» = A / 6 t * , . 

as determined by comparison of Eq. (17) and (15). 

For the general case, we have 

a 0 = R < 0 ) , a, = 1/2RM0), a a = -1/6R" (0) . 

One can expect a linear dependence of SV ( t ) for small 

when 


or 


a a t * / a i t » t / 2 1 0 < < 1 , 
t < < t. o / 3 „ 



~t ( m*) 


( 16 ) 

ahold 

(17) 

t h e 
(IS) 

( 1 9 ) 
i rn 0 ?!> 

;■ 0 ) 


Radiation-induced shift in threshold potential vs 


t :i. me .. 
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T h e f i r s t t h r e e p o i n t s s h owe d t h e 1 i n e a r d e p e n dence of S ' V ( t ) / d t <r> +• 
(30 min.) in both cases,, but the fourth point deviated from that 
dependence (Fig. 5). From this, we could conclude at once that 
t >3 1 = 30 min. In other word's, that t a was between one and two 
hours. We used the linear part of the curves (the first 2-3 
points) to find a and a,., and we then used the next part of the 
curve to find a». The constants a* and a 2 allowed us to calculate 
t 0 and A: 


■I/O a i / c 


a 1 ■ 


21 ) 


For the p-channel transistor (Fig. 6), we found t 0 to be 
approximately 1.1.0 min. and for the n -channel, t 0 70 min. 

F o r t. hi e p - c: h a n n el, t h e t h e o r e t i cal c u r v e , 

cVV (t ) /D t ot = a D -i- a,. t + a 2 t 2 , 

deviates from the experimental points only after 70 min., which is 
0 „ 6 4 1 , a n d for t h e n -channe 1 , t hi e d e v i a t i o n t a k e s p 1 ace a f t e r 4 5 
m i n . , w h i c h i s also 0 . 6 4 1 o . 



Fig. 6. Shift per unit total dose for p-channel transistor. 

For the n-channel (Fig. 7), we then plotted a more precise 
theoretical curve, adding one more term, a 3 t 3 : 

SV (t) /D^ot = a <:> + ait + a 2 t 2 + a 3 t 3 . 

The value of a 3 was found from Eq . (15) and (15a): 

a 3 = 3/2 a 2 2 /a = 0.00196 - 0.002. (22) 

We plotted this improved curve to t = 1 hr, which gives t/t D = 
0.86. There is no point in continuing the curve further, since 
the condition t/t D il must hold and the expansion (Eq. (17) is 



13 


ORIGINAL PAGE ?3 
OF POOR QUALITY 


d e f i n 1 1 e 1 y n o t v a 1 i d b e y o n d this point. Our 1 purpose here i s t o 
a n a 1 y 2 e t h e r e g i a n t /' t < 1 „ 



t( 


m/n 


) 


Fig. 7. Shift per unit total dose for n-channel transistor. 
Pure Annealing (t<<t 0 ! 

The case of pure annealing can be treated in the same way. 
Instead of Eq. (3) , we will use our general equation. 

Substituting Eq. (12) in Eq. (2) , we obtain 

SV = DtotER!0) + R’ (0) (t - l/2t r .) + 

+ 1/2R" <0> (t* - tt_ - l/3t r - 2 ) + 

+1/3R’ " (0) <t 3 - Zt^tr- + ttr- 55 - 1 /4tr- 3 ) 1 (23) 


or 


<5'V ( t ) /Dtot = R ( 0 ) - 1 /2R’ (0) t,_ - 1/6R" (0) t,- 35 - 1 / 12R 1 "' (0> t,. 3 + 

ER’ (0) - 1/2R" <0)t_ + R ? " (0)tr- 2 3t + 

( 1/2RR" (0) - 1/2R’ " (0) t,-) t 2 + 

1/3R’ " (0) t 3 . (24) 

Since we are concerned with the case t r -/t c ,<< 1, Eq. (21) can be 
s i mp 1 i f i ed : 

jV(t)/D tot = R (0) + R’(0)t + 1 /2R" (0) t 3 + 1 /3R ? " (0) t. 3 . (25) 

In particular, for the ln(t) dependence, we substitute Eq. (12) into 
E q . ( 22 ) an d ob t a i n , 

S"'J ( t. ) /DfccD-t = -C + A (t/t 0 ) -l/2A(t/t e3 )= + 2/3 A (t/t„) 3 (21) 

Then, as in the case of simultaneous irradiation and annealing, we 
can, in principle, determine a 0 , a x , a 2 from experimental data. 
We can then calculate t D and A: 

t. ca ~ ~l/2 a j. / a 2 


:L4 




A= 


< 28 ) 


Temperature and Dose-Rate Dependence of Response Function. 


It is worth noting that our lack of data on the dependence? of R<t) 
on dose rate and temperature at this point prevents development of 
the microscopic (quantum) theory of annealing. 

If we find that, in some cases, R(t) does not possess a ln(t.) 
dependence (Eq. (1)), our method still allows us to find the value 
of R(t) and its derivatives at. t = 0 through Eq. (5). Again, by 
conducting experiments for different dose rates (or temperatures) , 
one can obtain the dependence of RIO), R ' (0) and E" " (0) on these 
variables. Since R” (0) is proportional to l/t D and R" (0) is 
proportional to 1 /t Q 2 , we can determine the dependence of t„ on 
temperature or dose rate. 


III. LONG-TERM ANNEALING AFTER SEVERAL DISCRETE DOSES. 


A study of variable dose rate would gain for us more insight into 
the problem of annealing. Of course, it must be preceded by a 
detailed study of temperature and dose rate dependence of A, C and 
the characteristic time, t 0 „ 

We shall treat the specific: case of variable dose rate in a 
stepwise varying dose rate in detail elsewhere. Here we will 
treat the case of several impulse doses which have different dose 
rates (Fig. 6). By impulse dose, we mean doses whose duration, 1 t , 
is much less than the annealing time, i.e., 1 t<<t. We suggest 
that the response function has a ln(t> dependence, but each dose 
is described by its own character! st i c time, *t ra . We shall assume 
that. A and C are independent of dose rate. 

Thus, the integral response is the summation of the responses for 
e a c h d o s e r a t e s 


tr 

n 

,5 V ( t ) = Z J R ( II t -t =• 3 / 1 1 c ) Di dt ' ( 29 ) 

i 0 


For the long-term annealing in this case, we obtain 

ffV(t) = ERft/HoJDi, (29a) 

„ i. 

where D* = D i i t is the individual total dose for a dose rate D* 
and 1 t is the duration of the i-th impulse. 
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Fig. 8. Three discrete doses with different dose rates. 

For the ln(t/t i:> ) dependence, we obtain 

<5'V ( t ) = EC-C + A In (t/t D ) ]D ± . (30) 

i 

It is easy to show that Eq. (30) takes the form of a single impulse 
with the total dose: 


and with some effective character i sti c time, t Q . To prove? this 
statement and to find this effective char ac ter i st. i c time, it is 
sufficient to treat the case of two impulse doses. In this case, 
we have 


D 1 1 n ( t / 1 1 ) + D a ln(t/t a ) = D tcJt ln (t /tj . (32) 

Then Eq. (31) takes the form: 

(D,./D t ot ) (Da/D ■fc otr. ) 

lnC(t/ti) (t/t 2 ) "I = ln(t/t„), 


•from which 

(Di/D ) (D a /Dtot! 

to = (t i ) (t 2 ) . (33) 

Eq. (33) can obviously be generalised for any number of impulses: 

(Di/Dtot) 

to = TT ( 1 t o ) . (34) 

j. 

The principal conclusion to be drawn from E. (29) is that the 
contribution of a single "strong" k-th impulse, i . e. , an impulse 
whose dose rate is much larger than that of any other i—th 
impulse, does not produce much change in t ra , provided that its 
contribution to the total dose is small (D k <<D tot ) . 
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IV. NON-LINEAR EFFECT. 


We have assumed that the response is linear with an impulse dose 
and that the response -function depends only on the time lapse 
between irradiation and measurement. At this point, we will show 
h o w e :■! i s ting e x p e r i rn e n t a 1 d a t a s h o u 1 d b e e x a m i n e d to v e r i f y t h e s e 
t wo assump t i on s . 

We shall now examine nonlinear effects. The output,, that is, the 
shift. of the threshold potential,, is no longer assumed to be a 
1 i n e a r f u n c t i a n o f a n i m p u 1 s e dose, b u t. i s p r o p o s e d t. o b e 
proportional to an unknown impulse function, <5'F(D): 

cVV ( t ) = F: ( t - 1 •’ ) SF ( D ) . ( 35 ) 

The total shift, is then given by 


V* o fe = J R (t-t ) F M D ) d D , < 3 6 ) 

where F’ (D) - dF(D)/dD. In the? limit, in which t>>tr-, we obtain 

V ( t. , D tot ) = R < t ) F (Dtot ) - < 37 ) 

In contrast to the linear case, we now have not. one, but two 
u n k n o w n f unctio n s , R a n d F , but only one q u a n t i t y which c: a n b e 
experimentally determined: V. These two functions, however, 

depend on two different variables, t and D to t. This allows us, by 
fixing one variable, to determine the function which depends on 
t. h e at h er var i ab 1 e . 


V. 


T w o - V a r i a b 1 e R e s p a n s e Func t i a n , 


Next, we shall treat the case of a two-variable response function, 
R(t, t : ' ) . When t>>tr-, we obtain 


w 


tot R vt , O ) Dtot . 


(33) 


Eq. (12) gives F: only as a function of the first variable. We have 
devised a method to determine R as a function of the second 
variable. One needs to provide two impulse doses, D cl> and D <2> , 


separated by an interval t ’ 


(Fi 


9 ) 


W e t. hen have t h & f o 1 .1 o wing 


expression for the shift of the threshold potential: 

V ,a, (t) = R(t, 0)0 (1> + R(t, tOD' 2 ’, 


(39) 


provi ded t >>t „ 

Substituting the first term on the right-hand side of Eq. (13) 
Eq . (12), we ab t a i n 


1 > 

V tot 


< ZZ > 


( t ) 


( Zi > 


by 


R ( t , t =• ) 


Vt=,t < 1 > (t) /Dtot 


*! 


(40) 
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w h e r e V < 1 > a n d V (2> are t h (3 s h i f t s o f t h e t hre s h o 1 d p o t e n t i a 1 s 
res p e c t i v e 1 y -for t. he -first a n d s e c o n d impu .1 s e s 

S i n ce t i s a rbi t r ary, p r o v i d e d t > > t ’ , E q . < 4 0 ) d e t e r m i n e s R a s a 
•function of the second variable, Eq.(40>, together with Eq. (38),, 
compl etel v deter flu nes R as a function of two variables. 


V 



Fig. 9. Two discrete impulse doses. 


VI . CONCLUSION 

The phenomenological linear response theory developed here is 
capable of describing the entire range of the reversible kinetic 
process of irradiation/annealing. Both pure annealing and 
simultaneous irradiation/annealing have been treated in various 
time ranges. These include; small time intervals (t<<t 0 > , short- 
term annealing <t on the order of t,~) and long-term annealing 
(t >>t r -, t „ ) , as well as the saturation region ( <5'V approaching O) . 

A m e t h o d was also developed to determ i n e the c h a r a c t e r i s t. i c 
parameters of the response function. An important advantage of 
this method is that it is useful to determine the temperature and 
dose-rate dependence of these parameters. 

In dealing with reversible linear kinetic process, one is 
naturally interested in the limits of reversibility. To study 
these limits and to examine the transition are between reversible 
and irreversible processes, we have presented the case of two 
variable response functions and have treated a response which is 
non-linear with respect to an impulse dose. 

It was pointed out that in the case of simultaneous irradiation 
and an n e a 1 i n g , e specially a t. t he limit s o f s m a II time, t< <. 1 0 , * n d 
large time, t >>t 0 , the respo n s e per u n i t total d ose. A ' V ( t ) / D t . <=>*, 
is more convenient to deal with than the shift itself. 

Finally, we treated the case of several different, impulse doses 
and showed that they act as a single impulse dose with a common 
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0 f f e c t i v 0 c: h a r a c t e r i s t i c t i m e - I n p a r t i c. u 1 a r , i t w a s s h o w n t h a t a 
given impulse dose, if it makes a negligible contribution to the 
t. a t. a 1 do s e , does not change the eft e c t i v s c h a r a cter i s t i c: t i m e . 

S u c h a c o n c: 1 u s i o n i s v a 1 i d , p r o v i d e d t h a t t h e i m p u 1 s e d o s e d o e s 
n o t: e x c e e d t h e 1 i m i t s o f i r r e v e r s i b i 1 i t. y f o r t h e w h o 1 e p r o c e s s „ 

We must realise that our main goal is still to construct a 
microscopic theory,, Only' quantum theory can give us an insight 
into the complicated processes related to radiation damage and 
a n n e a 1 i n g a n d p r e d i c: t tern p e r a tur e a n d d o s e r a t e d e p e n d e n c e . 

To build such a theory, however, will be a difficult task, despite 
the fact that solid-state quantum theory has reached quite a 
s a p h i s t i c a t e d s t a t e . T h e m a i n r e a s o n is t h a t h e r e w e a r e 
c o n c e r n e d w i t h k i n e tic p r o c e s s e s , w h e reas t h e p r e s e n t q u a n t. u m 
t. h e o r y stu d i e s m o s 1 1 y equi 1 i b r i u m o r qu a s i - e q u i 1 i br i u m pro c e s s e s . 
In addition, we are concerned mainly' with processes far from 
equilibrium, and even with the transition to irreversible 
processes , 

To make any progress toward a foundation of this microscopic 
theory, we need, as a necessary first step, to have a general 
phenomenological theory.. The linear response theory outlined in 
this paper can serve as an important, step toward a quantum theory 
of r a d 1. a t i o n d a m a q e a n d a n n e a 1 in g „ 

I n c o n c 1 u s i o n , w e ha v e 1 a i d t h e g r o u n d w o r k f or a n a 1 y t i c a 1 .1 y 
characterizing the radiation and annealing response of devices of 
both variable dose and dose rate, which will greatly aid in the 
s I m u 1 a t i o n o f s p a c e r a d i a t i o n e f f e c t s . 


V 1 1 „ RECDMMENDAT I ONS 

We regard as the first priority the study of the temperature 
dependence of the characteristic time, t D , and the parameters A 
and C. 


2„ Additional directions for experimental study' are: 
a. Dose-rate dependence 


b ., The 

sat ura t 

ion region, in 


w h i c h 

anneal i ng 

approaches 

comp 1 et i on 

( V app 

r o a c h ee s z e r a ) a n d 






c „ The ir 

radi at i 

o n c: o n d :i. t i a n s b e y 

on 

d whicl 

I t h e k i 

n et i 

c processes 

become irre 

versi b 1 

e and the extent 

of 

these 

cond i t i 

ons . 
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